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ABSTRACT 

This paper aims to consider the evidence that debris disks are self-stirred by the 
formation of Pluto-size objects. A semi-analytical model for the dust produced 
during sclf-stirring is developed and applied to the statistics for A stars. We 
show that there is no significant statistical difference between fractional excesses 
of A-stars <50Myr old, and therefore focus on reproducing the broad trends, 
the "rise and fall" of the fraction of stars with excesses that the pre-stirred 
model of Wyatt et al. (2007b) does not predict. Using a population model, we 
find that the statistics and trends can be reproduced with a self-stirring model 
of planctcsimal belts with radius distribution N(r) cx r -0 8 between 15-120 AU, 
with width dr = r/2. Disks must have this 15 AU minimum radius in order to 
show a peak in disk fraction, rather than a monotonic decline. However, the 
marginal significance of the peak in the observations means that models with 
smaller minimum radii also formally fit the data. Populations of extended disks 
with fixed inner and/or outer radii fail to fit the statistics, due mainly to the 
slow 70 fini evolution as stirring moves further out in the disk. This conclusion, 
that debris disks are narrow belts rather than extended disks, is independent of 
the significance of 24 /im trends for young A-stars. Although the rise and fall 
is naturally explained by self-stirring, we show that the statistics can also be 
reproduced with a model in which disks are stirred by secular perturbations 
from a nearby eccentric planet. Detailed imaging, which can reveal warps, sharp 
edges, and offsets in individual systems is the best way to characterise the stirring 
mechanism. From a more detailed look at (3 Pictoris Moving Group and TW 
Hydrae Association A-stars we find that the disk around (3 Pictoris is likely the 
result of secular stirring by the proposed planet at ~10 AU; the structure of the 
HR 4796A disk also points to sculpting by a planet. The two other stars with 
disks, HR 7012 and r\ Tel, possess transient hot dust, though the outer r\ Tel disk 
is consistent with a self-stirred origin. We suggest that planet formation provides 
a natural explanation for the belt-like nature of debris disks, with inner regions 
cleared by planets that may also stir the disk, and the outer edges set by where 
planetesimals can form. 

Key words: circumstellar matter - stars: planetary systems: formation - stars: 
planetary systems: protoplanetary discs. 



1 INTRODUCTION 

Debris disks are the disks of dust found around nearby 
main sequence stars through their thermal emission (for 
a recent review see Wyatt 2008). The dust itself is short- 
lived compared to the lifetime of the star so is be- 
lieved to be continually replenished through collisions 
between km-sized planetesimals (Wyatt & Dent 2002; 
Quillen et al. 2007), much in the same way that dust 
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in the zodiacal cloud is replenished through collisions in 
the asteroid belt and Kuiper belt (Dermott et al. 2001; 
Moro-Martm & Malhotra 2003). While there remain dif- 
ficulties in growing dust from its initially sub-/im size 
into >km-sized planetesimals (Blum & Wurm 2008), the 
widely invoked coagulation and core accretion models for 
planet formation rely on relative velocities in the proto- 
planetary disk being low enough that collisions result in 
net accretion, rather than destruction. As it seems that 
the opposite is the case in debris disks, these disks must 
have been stirred at some point so that collisions occur at 
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>l-10m/s, which typically corresponds to eccentricities 
and inclinations for the disk material of e > 1CP 3 to 10 -2 
(e.g. Kcnyon & Bromley 2008). 

To explain debris disks it is not normally necessary 
to understand how (or when) they were stirred. It is suffi- 
cient to invoke a "pre-stirred" debris disk — one that was 
stirred when the star was born. For example, it is possible 
to explain the statistics of dust found around A-stars of 
ages >10Myr by assuming that all stars are born with a 
pre-stirred planetesimal belt that evolves due to steady- 
state collisional erosion (Wyatt et al. 2007b). The diver- 
sity seen at different ages then reflects their different ini- 
tial masses and radii. A similar conclusion was reached 
for debris disks around sun-like stars (Lohne et al. 2008). 

However, recent results on the presence of debris 
disks around young A-stars are challenging this view by 
showing that the fractional excess at 24 (jm from hot dust 
increases from 3 Myr (the time at which most protoplane- 
tary disks have dissipated) to a peak at ~ 10-30 Myr, fol- 
lowed by a slow decline as expected by steady-state evolu- 
tion (Hernandez et al. 2006; Currie et al. 2008a,b). This 
peak has been interpreted as evidence for self-stirring, 
where debris are created when the largest objects become 
massive enough to stir planetesimals to fragmentation ve- 
locities. 

Self-stirring models follow the evolution of an ex- 
tended planetesimal belt from the protoplanetary disk 
stage, allowing both the size distribution and the velocity 
distribution to evolve in a self-consistent manner. These 
models find that stirring of the planetesimal belts occurs 
when planets reach Pluto-size, which depends strongly on 
their distance from the star, as well as on the mass surface 
density of solid material there: tpiuto oc P/E, where P is 
the orbital period and E is the surface density of plan- 
etesimals (e.g. Lissauer 1987). For a typical disk model 
E = Eo r~ s , so tpiuto depends strongly on radius r: 

tpiuto oc r 3 / (Eo y/M^j . (1) 

Thus, compared to pre-stirred disks, self-stirring means 
that farther disk regions are stirred at later times because 
it takes longer for Pluto-size objects to form there. This 
evolution means that individual disks can show a peak in 
emission at ~10-30Myr, set by the time to form Pluto- 
size objects at the inner disk edge (Hernandez et al. 2006; 
Currie et al. 2008a; Kenyon & Bromley 2008; Wyatt 
2008). That is, the observed evolution occurs if debris 
disks tend to have ~10 AU inner holes. The overall dust 
content at <10 Myr is low because such times are needed 
for Plutos to form at the inner edge of the disk. This lack 
of dust implies that accretion with minimal fragmenta- 
tion is ongoing in young (<10Myr) debris disks and the 
outer regions of older debris disks. 

Though the Kenyon & Bromley (2008) self-stirring 
models are qualitatively consistent with the ^10-30 Myr 
peak in fractional excesses, there has been no quantitative 
test of whether these models can reproduce the observed 
A-star statistics, or how the observations constrain disk 

1 Self-stirring is also called delayed stirring in the literature, 
but following Wyatt (2008) we consider that what we are call- 
ing self-stirring is a sub-set of possible delaycd-stirring models. 



parameters. In this paper, our aim is to address this is- 
sue using the Wyatt et al. (2007a) debris disk evolution 
model, modified to include self-stirring (as outlined in 
Wyatt 2008). 

In particular, there are several issues to address. Be- 
cause the peak in excesses occurs when Plutos form at 
the inner disk edge, the timing of this peak is highly de- 
pendent on the radius of the cleared region and on mass 
surface density (Eq. 1). Yet we know that protoplane- 
tary disks have a range of masses and surface densities 
(e.g. Natta et al. 2000; Andrews & Williams 2005, 2007), 
and it seems unlikely that the inner clearing would be 
at the same radius for all disks. Therefore, we wish to 
find whether self-stirring models can reproduce the A- 
star statistics, and if so, whether they put a strong con- 
straint on disk inner radii as Equation (1) suggests. 

Another issue is the extent of debris disks. 
Kalas et al. (2006) note that disks resolved in scattered 
light appear to be either relatively narrow belts 20-30 AU 
wide, or extended disks wider than 50 AU (however, a 
disk that appears extended may result from blowout 
of small grains created in a relatively narrow planetes- 
imal belt). In fitting the statistics for (>10Myr) A-stars 
Wyatt et al. (2007b) considered disks to be the former; 
belts centred at radial distance r with an assumed width 
of dr = r/2. However, self-stirring allows the interesting 
possibility that disks similar in extent to observed proto- 
planetary disks — from near the star to many hundreds 
of AU — may evolve to look like narrow belts because 
only regions of recent Pluto formation may be luminous 
enough to be detected. Therefore, we also wish to evalu- 
ate whether debris disks tend to be "narrow belts" or "ex- 
tended disks" (we use these terms consistently through- 
out to refer to these types of disks). 

A final issue is whether debris disks could be stirred 
by something other than the formation of Pluto-sized ob- 
jects. Though self-stirring has been the only proposed 
mechanism, secular perturbations from planets (another 
subset of delayed stirring) seems to be an equally vi- 
able way to stir debris disks (Mustill & Wyatt 2009). 
This hypothesis is partly motivated by stars with debris 
disks known or predicted to harbour planets, such as Fo- 
malhaut and /3 Pic (Mouillet et al. 1997; Quillen 2006; 
Kalas et al. 2008; Lagrange et al. 2009a). The question 
is therefore whether we can identify the more important 
stirring mechanism, both at a population level and for 
individual objects. 

The layout of this paper is as follows. In §2 we outline 
the Wyatt et al. (2007a, b) model, and the modifications 
made to include self-stirring. We empirically fit some 
model parameters to reproduce the self-stirring models of 
Kenyon & Bromley (2008) in §3. As self-stirring results 
in not only excess evolution, but also spatial evolution as 
Plutos form farther from the central star, we show how 
disk surface density profiles evolve and vary with model 
parameters. In §4 a population model is compared with 
the statistics for A-stars. We briefly consider a planet- 
stirred population model, take a more detailed look at 
resolved /3 Pictoris Moving Group and TW Hydrae A- 
stars in the context of delayed stirring, and discuss other 
influences on debris disk structure in §5. We summarise 
our main conclusions in §6. 
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2 ANALYTICAL SELF-STIRRING MODEL 

This section describes the analytical model. It is an ex- 
tension of the model described in Wyatt et al. (2007a, b) 
to include delayed stirring by splitting disks into a series 
of concentric annuli. The implementation uses the equa- 
tions in §2 of Wyatt et al. (2007a, b) and sums over 100 
logarithmically spaced annuli to create a disk with radial 
extent. The self-stirring prescription is semi-empirical, 
based on the more detailed models of Kenyon & Bromley 
(2008) (see §3.1). Here, we briefly summarise the model 
and refer the reader to §2 of Wyatt et al. (2007a, b) for 
details omitted here. 

We use a Minimum Mass Solar Nebula (MMSN, 
Weidenschilling 1977) surface density profile to specify 
the disk mass 



a tot = 3.5 x 10- L 'K(3q - 5)" 1 (HT 9 D b i) 5 ~ 39 



(3) 



E = r\ M* E r~ 



(2) 



where r\ is a scaling parameter reflecting a range of disk 
masses, h'h is in Solar units, radial distance r is in 
AU, and the surface density at 1AU is Eo = 30 gem -2 
(1.1 /AU 2 in our units). The power-law index 8 is typ- 
ically 1—1.5; we use 1.5 as our canonical value. This sur- 
face density provides roughly the amount of solid material 
contained in the outer Solar System. The surface density 
scales linearly with stellar mass, consistent with mm ob- 
servations of protoplanetary disks. However, the range 
of stellar masses in our model is much smaller than the 
expected range of 77, so the factor M* is relatively unim- 
portant (Natta et al. 2000; Andrews & Williams 2005). 
For extended disks the inner and outer radii are inde- 
pendent and specified by n n and r ou t- Narrow belts are 
specified by their radii r m id and width dr (so for belts 
nn,out = r mi d ± dr/2). 

Where Wyatt et al. (2007b) considered a range of 
total masses M tot , we consider a log- normal distribution 
about rjmid- 2 There is as yet no evidence of a positive 
correlation between disk mass and radius (which might 
argue for the distribution of 77 to be narrower than that 
for Mtot), and thus we use the la width of 1.14 dex ob- 
served for disk masses in Taurus (Andrews & Williams 
2005). 

The planetesimal disk is assumed to be in colli- 
sional equilibrium with a size distribution in each annu- 
lus defined by N(D) = KD 2 - 3q , where K is a constant, 
q = 11/6 for an infinite collisional cascade (Dohnanyi 
1969), and D is planetesimal diameter. That distribution 
is assumed to hold from the largest planetesimal in the 
disk, of diameter D c (in km), down to the size below 
which particles are blown out by radiation pressure as 
soon as they are created, Dm (in /an). If q is in the range 
5/3 to 2 then most of the mass is in the largest planetes- 
imals while the cross-sectional area atot is dominated by 
the smallest particles: 



2 Previously each disk was assigned a total mass distributed 
about M m id independent of disk location. This led to typi- 
cal disks at r = 120 AU having about 6 times lower So than 
disks at 3 AU (for <5 = 3/2). However, this difference is smaller 
than the dispersion in surface density and disk mass so is not 
particularly important. 



in AU 2 . 

We assume that particles act like blackbodies, so the 
fractional luminosity of grains / = Lir/L 4 = a to t/47rr 2 
for the annulus at r. The grain blowout size is 



D b i = 0.8 (L*/Af*) (2700/p) 



(4) 



where Db\ is in /im, L* and M* are in Solar units, and 



density p is in kgm 
worked out from: 



T 



278.3 L°- 25 r- 05 



The dust temperature can be 



(5) 



We also assume that the central star acts like a black- 
body. 

Because we use a fixed size distribution with q = 
11/6 in our model, the long-term evolution of the disk 
is determined by the removal of mass from the top end 
of the cascade. Planetesimals have a disruption thresh- 
old Qn, and eccentricity e. The collisional lifetime of the 
largest planetesimals of size D c at a radius r in an MMSN 
disk is 

i = 2.2 x l - 10 r 23 / 6 D c Q* D 5/6 e" 5/3 M~ 7//3 Eq 1 77- 1 (6) 

in Myr, obtained by substituting E = M to t/2Trrdr into 
Wyatt et al. (2007b) Equation (9). We have simplified 
the relation here by assuming that X c = D cc /D c -C 1, 
where D cc is the smallest planetesimal that has enough 
energy to catastrophically destroy a planetesimal of size 
D c . This condition applies for the e > 0.01 eccentricities 
found for self-stirring models (Wyatt et al. 2007a). 

Assuming that collisions are the only process affect- 
ing the evolution of the surface density, then 

S (rt)-i E(r '° } t<tstil (7) 

Mr,i; \ E(r,0)/[l + (*-t stir (r))/t c (r-,0)] t > i Btir [<> 

where E(0) is the initial surface density profile, t s tii is 
the delay until stirring (assumed to be in Wyatt et al. 
2007b), and t c (r, 0) is the collisional lifetime at r at that 
initial epoch. Where t < t a ti r , the initial surface density 
profile applies. 

Because more massive disks process their mass faster 
(t c oc I/t?), the surface density at late times (t 3> t c and 
t > tstir) a t a given radius is independent of initial surface 
density. Written in terms of the surface density of cross- 
sectional area in the disk (or simply "surface density" 
because r e fl oc E, T e g = a to t /2irrdr — 2fr/dr). 



T e f£,i 



1.2 x 10" 9 r 7/3 D° A VD 
x e-^M-^L-^it-t^y 1 



- 5/6 



(8) 



This prescription means that for a disk of known size 
at any given age, there is a maximum infrared luminos- 
ity, / ma x, that can remain due to collisional processing. 
While / ma x can be applied to unresolved stars based on 
the 24 — 70 /im colour and some assumption about disk 
extent (Wyatt et al. 2007a), r e ff, m ax necessarily requires 
resolved observations for comparison. We return to the 
application of r e ff, m ax to self-stirred disks in §3.3. 

The final component of the model is the implementa- 
tion of self-stirring, where successive annuli can be stirred 
at later times with increasing r (such as Eq. 1). Other 
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mechanisms for delayed-stirring, such as secular pertur- 
bations from an eccentric planet (e.g. Mustill & Wyatt 
2009), can be implemented with different t 3 tiv (see §5.1). 
The introduction of self-stirring also requires us to spec- 
ify the level of disk emission before a disk is stirred. We 
assume that the emission is proportional to the surface 
density, but reduced by a factor Xdciay with respect to 
that expected from Equation (3) because the collisional 
cascade has not yet begun. We estimate Zdday through 
comparison with the Kenyon & Bromley (2008) models 
below. 



3 DEBRIS DISK EVOLUTION FOR 
INDIVIDUAL STARS 

In this section we consider model evolution for in- 
dividual disks, first making comparisons with the 
Kenyon & Bromley (2008) models, and then illustrating 
the surface density evolution. We then consider the impli- 
cations for transience in the context of the Wyatt et al. 
(2007a) model and the effects of disk extent on excess 
evolution. 

3.1 Comparison with Kenyon & Bromley 
models 

We compare our model with the Kenyon & Bromley sim- 
ulations of planet formation and debris disk evolution 
around A-stars to derive an empirical self-stirring time 
and idoiay that reproduces their excess evolution. Though 
Kenyon & Bromley (2008) derive relations for the time to 
reach the peak excess (their Eq. 55), and the radial time 
dependence for forming 1000 km objects (their Eq. 41), 
they do not derive a relation for the peak excess at differ- 
ent radial locations. Because our aim is to reproduce their 
excess evolution, we use the following empirical relation. 

t„tir = 4 x 10 -4 r 3 tT 1/2 M~ 3/2 (9) 

in Myr. We first derived the numerical pre-factor and the 
exponent of r\ in this relation with a \ 2 minimisation 
procedure over models with 1.5 < M* < 3 and 1/3 < 
rj < 3 at 24 and 70 (im. We fix the radial and stellar 
mass dependence because these are set by the orbital 
period and our assumption of E a M*. However, this 
formal method ignores some important model differences, 
in particular the continued stirring that appears to cause 
the Kenyon & Bromley (2008) models to decay faster at 
70 /zm (discussed below). We therefore arrive at Equation 
(9) by focusing on the 24 fim peak excesses. This different 
emphasis changes the values obtained in the formal fit at 
24 fxm by ~10-20%. 3 

Though the dependence on rj in Equation (9) is 
weaker than implied by Equation (1) (—1/2 vs. — 1), 

3 Similar variation is also found depending on how the 
Kenyon & Bromley models are weighted in the minimisation 
procedure. For example, lower weights for smaller excesses 
makes the peak excesses fit better, whereas constant weights 
fit the late time evolution better but underestimates the peak 
excesses. Allowing £dclay to vary results in values in the range 
~0.01-0.02. 



the two should not necessarily agree. Kenyon & Bromley 
(2008) find that the time to peak dust production is 
oc ry -1 M^ -1 ' 5 (as predicted by Eq. 1), but that the time to 
peak luminosity is oc r)~ 2 ^ Af" 1 and takes much longer 
(their Eqs. 48 & 55). In comparing these proportional- 
ities, it is important to remember that the peak dust 
production and peak luminosity are linked, but not di- 
rectly. The observed luminosity depends on the history 
of the dust production rate and the rate at which small 
grains are being removed by radiation forces. Further, 
our empirical relation sets the evolution at specific r, 
rather than describing the global evolution of specific disk 
properties. Though Equation (9) is derived largely from 
the 24 /im evolution, we reproduce the / evolution, from 
which the Kenyon & Bromley peak luminosity is derived 
equally well. Our scaling of t s y r with r\ is probably weaker 
than their —2/3 because we set the onset of stirring, not 
the peak excess and higher mass disks reach their peak 
excess more rapidly once stirring starts (Fig. 1). We refer 
to the radial location where stirring has just begun (i.e. 
where t = t st ir) as r s tir- 

Figure 1 shows a comparison of the 
Kenyon & Bromley model for an A2 (2.5 Mq) star 
with an equivalent one from our model. The 24 /im ex- 
cess increases to a peak around 10 Myr, and then declines 
(see Hernandez et al. 2006; Currie et al. 2008a; Wyatt 
2008; Kenyon & Bromley 2008, for comparison with 
observations). In this example, the disk extends from 
30-150 AU, and has r? = 1/3, 1, and 3. We set the level of 
disk emission before stirring to a^doiay = 0.01, which gives 
a reasonable match to the excesses prior to the formation 
of Pluto size objects. We take other model parameters 
from the best fit model of Wyatt et al. (2007b); e = 0.05, 
D c — 60 km, and Qfj = 150 J kg . Though the largest 
objects in the disk are roughly Pluto size when the 
collisional cascade begins, Kenyon & Bromley find that 
objects larger than ~ 1-10 km continue growing. Thus 
the effective maximum planetesimal size D c for our 
model is of order 10 km because the smallest dust derives 
from a reservoir of objects of this size. Different choices 
for D c , e, and Qd change the rate at which excesses 
increase after the disk is stirred, though still give peak 
excesses that agree with the Kenyon & Bromley models 
within factors of a few. 

To achieve the agreement shown in Figure 1 we 
must temporarily modify our model to match some 
Kenyon & Bromley assumptions about dust properties; 
they set Dbi = 1 /im, and use a "greybody" emission law 
Fy oc Bi/(T) (l — exp(A /A)) with A = 10/xm. They use 
Yi et al. (2001) luminosities, somewhat brighter than the 
Schmidt-Kaler (1982) luminosities used in our model. Fi- 
nally, to ensure a good match between models we apply 
a factor 2.5 decrease to the Kenyon & Bromley excesses. 
The need for this factor is unsurprising considering the 
model differences and could, for example, arise from dif- 
ferences of a few percent in the average size distribution 
index q between D c and Dbi- With these differences taken 
into account, both models produce fairly similar results. 
The delay in the rise of 24 and 70 /im excesses is similar, 
as are their magnitudes over the range A0-F2 (~1.5- 
3M ) and rj = 1/3-3. 
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Figure 1. Evolution of 24 /an (left) and 70 fim (right) excesses for r\ = 1/3, 1, and 3 for a disk from 30 to 150AU from our model 
(solid lines) compared to models from Kcnyon & Bromley (2008) (dashed lines), for an A2 (2.5 Mg) star. The main difference 
between the models at 70 ^m is caused by continued accretion and subsequent stirring in the Kenyon & Bromley (2008) models 
(see text). Also shown is the effect of the increasing main-sequence luminosity on an r\ = 1 disk with Dbi = 1 A* 111 (s re y dotted line) 
and then also allowing £>bi to vary (grey dot-dashed line). The overall effect of stellar evolution due to changing Z?bi and hotter 
grains (i.e. dot-dashed line) is a small increase in 24 fim excesses and a similar decrease in 70 /im excesses. 



The main difference is the decline in 70 /im excesses 
at late times. Our model shows a fairly slow decline until 
after 200-600 Myr, when stirring reaches the outer disk 
edge and the excess decreases more rapidly because mass 
is being lost at all radii. The Kenyon & Bromley (2008) 
models show a faster decay. In their model, the largest 
objects continue to accrete once they reach Pluto size, 
thus increasing the rate of decay through removal of mass 
from the collisional cascade. More important is the in- 
creased stirring due to their continued growth, and the 
subsequent higher collision rates (S. Kenyon, priv. com- 
munication). This importance can be understood from 
the strong e dependence in Equation (6). 

Figure 1 also shows the effect of including stellar evo- 
lution in our model (scaled to show relative differences). 
For fixed Dm , increasing L* increases grain temperatures 
and the 24 /im excess is higher. If Dbi is allowed to vary, 
the disk emission drops slightly due to smaller a tot as 
Dbi increases and the overall effect is a small increase in 
24 fim emission and a similar decrease at 70 fim. 

After a few hundred Myr, the Kenyon & Bromley 
70 fim excesses also show a break towards faster de- 
cay. They attribute this decrease to stellar evolution as 
the star nears the end of its main-sequence lifetime (of 
650 Myr for 2.5 Mq). Higher stellar luminosity increases 
the importance of PR drag relative to collisions and the 
disk emission drops more rapidly as grains spiral toward 
the star. However, any mass lost to PR drag is probably 
due to the assumption of a fixed Dbi, as these models 
are of observable disks, which are generally not tenuous 
enough to suffer PR drag before grains collide (Wyatt 
2005). For the A2 star in Figure 1, Dm should increase 
from about 10 to 20 fim during the MS lifetime, so fix- 
ing Dbi = 1 fim allows PR drag to remove grains that 
would have instead been blown out of the system. Also, 
a disk is unlikely to become PR dominated by realistic 



increases in stellar luminosity. 4 Our model shows a break 
at 200 Myr for the rf — 3 (top) line at 70 fim when stir- 
ring reaches the outer edge of the disk. Given that PR 
drag is unlikely to be the cause, the further drop in the 
Kenyon & Bromley models after ~200Myr may in part 
be for the same reason. 

In summary, we have shown that our simplified 
model is in good agreement with the Kenyon & Bromley 
(2008) models, with the main remaining discrepancies 
arising for older stars at 70 fim. The effect of changing 
main-sequence luminosities is small, and continued accre- 
tion by the largest bodies leads to smaller 70 fim. excesses 
after the peak is reached. We retain our prescription for 
-Dbi, the main cause of the initial differences because it is 
well established that grains should respond to radiation 
forces (e.g. Burns et al. 1979). We allow x dc i a y to vary, 
though it is usually unimportant because most emission 
comes from stirred regions. 



3.2 Evolution of surface brightness 

Debris disks around some of the closest stars have been 
resolved, providing spatial information that cannot be 
derived from unresolved photometry. In this section we 
show disk profiles derived from our model for comparison. 
With delayed stirring there are two possible ways for the 
disk to evolve when stirring begins. These are set by the 
collisional (t c ) and stirring (i s tir) times. These two modes 



4 For fixed M+ and r the orbit decay timescale for smallest 
grains (which are most affected by PR drag) is r.pR <x Z?bi/£*. 
However, because Z?bi increases with L* (Eq. 4), the PR 
timescale for the smallest grains is constant. The collisional 
lifetime of the smallest grains increases slowly as t co n <x y/ Db\- 
For fixed f3 = 0.5, the relative importance of PR drag vs. col- 
lisions ?7o = tpa/teoll ok 1/vT*, and therefore changes little 
for the factor ~few changes in main-sequence L*. 
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Figure 2. Evolution of r cft for sclf-stirrcd (solid lines) and pre-stirrcd (dotted lines) disks at 1, 10, and 100 Myr. In the left panel 
tc 3> istir so evolution is "slow" and similar to the pre-stirred case. In the right panel t c -C t s t; r so the evolution at r st ; r is "fast" 
and much more violent as the disk reverts to the equilibrium state. Lines are offset slightly to remove ambiguities. The disks 
extend from 5— 150 AU around an A2 star and have parameters: rj = 1, = 150 J kg , e = 0.05, aidclay — 0.01, and D c = 60 km 
(left panel) and D c = 1km (right panel). The peak surface density is larger for the fast self-stirred disk because more mass is 
concentrated in smaller size objects. 



of evolution are shown in Figure 2, which shows the r e e 
evolution vs. radius in an extended disk at several times. 

In Figure 2, we also compare the evolution with a 
pre-stirred disk — one that is stirred when the star is born 
at t = (dotted lines). In this case the surface density 
profiles increase with radius to a broad peak. The slope of 
the inner region, where mass at all radii is being depleted 
collisionally, is set by r e ff, m ax <x r7//3 (Eq. 8). Outside the 
peak the primordial T e g remains because the collisional 
time is longer than the disk age there (i.e. r fr, ma x oc 
r -3 / 2 ). The location of the peak moves outward with time 
(rcxi 6 / 23 ,Eq. 6) as more distant regions begin to decay. 

The solid curves in Figure 2 (left panel) show the self- 
stirred evolution when the collisional time is longer than 
the stirring time. Compared to the pre-stirred case, the 
difference is the XdeUy drop outside r st i r where the disk 
has not been stirred and accretion is ongoing. For this 
disk little decay happens at r s ti r immediately, but begins 
later when t > t c . This evolution is therefore similar to 
the pre-stirred case, and we term this mode of evolution 
"slow" for a self-stirred disk. Compared to a pre-stirred 
disk, the factor :£dciay leads to a narrower observed annu- 
lus if the disk is too faint to be detected outside r st i r . For 
comparison, see Figure 9a of Kenyon & Bromley (2004), 
which shows a qualitatively similar profile. 

In the other limit, when the collisional time is shorter 
than the stirring time (Fig. 2, right panel), for a given age 
the disk at r st i r has much more mass than it would if it 
were pre-stirred (the same mass as at t = 0). The onset 
of stirring is therefore violent, with the rapid shedding of 
mass just inside r st i r as the disk reverts to its "equilib- 
rium" state, that of a pre-stirred disk where the decay is 
independent of initial mass. This evolution was hinted at 
by Dominik & Decin (2003), where the excesses of disks 
with longer delay times decayed faster with time (their 
Fig. 3). We call this evolution the "fast" mode. These 
disks appear transient because the surface density can 



be significantly above the expected pre-stirred level for 
their age (see §3.3). 

However, we consider fast self-stirred disks with very 
sharp surface density profiles unlikely, because disks with 
short collisional times will start to decay before Pluto-size 
objects form. Planetesimal eccentricities increase as the 
largest objects grow, and t c decreases accordingly. Thus, 
a disk that would stir in the fast mode in our model 
would in fact begin to decay earlier when t > t c (unless 
planetesimal eccentricities increase rapidly, in which case 
the time difference is small). 5 Because the largest objects 
continue to grow after the disk is stirred, eccentricities 
continue to increase. After stirring, t c continues to de- 
crease until the disk begins to decay. Thus, we also do 
not expect disks where the stirring time is significantly 
shorter than the collision time. 

By equating the stirring and collisional timescales 
(Eqs. 9 & 6), we can derive the conditions required for 
the boundary between these regimes of evolution at r a tir 
(assuming X c <C 1). Evolution will be "slow" at r st h- if 

U = r 5 J° D c Q* D 5/(i e- 5/3 M~ 5/6 ^ > 2 x TO 6 . (10) 

At 50 AU on the 10 Myr curves, the disks illustrated in 
Figure 2 have TZ = 4 x 10 6 (left panel) and 6 x 10 4 
(right panel), 2 times greater, and 33 times smaller than 
Condition (10) respectively. Because the collisional time 
increases more strongly with distance than the stirring 
time, inner disk regions are more prone to fast evolution, 
as are disks with weak or small planetesimals. There is 
some tendency towards fast evolution for more massive 
disks, though if the stirring time were simply the Pluto 
formation time (Eq. 1) this condition would be indepen- 
dent of surface density. Because disks with TZ ~ 10 6 are 
the most physically plausible within our model, we check 

5 In §5.1 we show that disks that evolve in the fast mode 
are possible if the stirring mechanism is secular perturbations 
from an eccentric planet. 



© 0000 RAS, MNRAS 000, 000-000 



Are debris disks self-stirred? 7 



the distribution of 1Z when generating population models 
in §4. 

The evolution of slow self-stirred disks makes clear 
predictions for resolved observations of debris disks. 
While the brightest disks will yield broad surface density 
profiles, fainter disks only have an annulus of detectable 
emission because the surface density drops away interior 
and exterior to where the peak emission occurs. The ob- 
served width of the disk depends on how fast the disk is 
decaying relative to its age; however, we do not expect 
surface density profiles as sharp as shown in the right 
panel of Figure 2 for self-stirred disks. Of course, the 
disk extent is also set by where planetesimals can form 
and where they are not affected by external effects such 
as clearing, accretion, or ejection by planets. 

For slow self-stirred disks the peak emission is set by 
the collisional time (i.e. the peak surface density in the 
left panel of Fig. 2 is slightly interior to r s _ir at 100 Myr) 



If the difference between the level of emission at the peak 
and at r st i r is small (as in the left panel of Fig. 2 at 1 Myr), 
and the disk beyond r st i r is too faint to detect, then the 
outer edge of slow self-stirred disks will still appear to 
increase with time as r oc t 1 ^ 3 (Eq. 9). 
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Figure 3. Evolution of narrow belts and extended disks (as 
labelled) in / vs. T2A— 70 space from lMyr to lGyr. Disk 
_!4— 70 increases as stirring moves outward (so evolution is 
left to right). Here, / is derived from a blackbody at r24-7_ 
(i.e. what detection at 24 and 70 fim would imply). Solid 
lines show 24 and 70 fim blackbody detection limits from 
Wyatt et al. (2007b). Extended disks have ri n _ out = 15- 
120 AU and if = 0.05 - 1.05. For belts r mid = 20, 30, 40, 
50, 60, and 70 AU, dr = r m ;_/2, and 77 = 1. Belts increase in 
radius from left to right, and extended disks increase in / with 
higher r\. Other model parameters are as in the left panel of 
Figure 2. 



3.3 Transience 

In the Wyatt et al. (2007a) model, narrow pre-stirred 
planetesimal belts have a maximum luminosity. This 
property arises because the decay time depends on the 
disk mass, and all disks were assumed to be stirred at 
t = 0. With an assumption about their width dr, the 
value /max can be calculated for unresolved belts with 
observations in several bands, which gives an estimate of 
their radii. If the disk luminosity is significantly higher 
than / ma x the excess is unlikely to arise from steady-state 
evolution. 

Ideally, comparisons would be made between re- 
solved disks and the predicted maximum surface density 
Toff.max- As discussed above, we do not expect self-stirred 
disks to have r e g significantly greater than r e ff ,ma X ; those 
that do may be transient. However, because disks may 
be secularly stirred in the fast mode (§5.1), they can ap- 
pear transient without being so in the sense meant by 
Wyatt et al. (2007a). Another issue, whether disks are 
extended or narrow belts is important because disks that 
are truly transient may not appear to be if they are as- 
sumed to be wider than they actually are. Therefore, if 
disks are narrow belts then /max is a reasonable indicator 
of transience. 



3.4 Radius evolution 

There are considerable differences between the evolution 
of narrow belts and extended disks. Figure 3 shows the / 
vs. r24-7o evolution for disks extending from 15-120 AU 
for a range of initial surface densities, and a series of 
narrow belts with the same surface density centred at 
different r m id between 20 and 70 AU (r24-7o is the radius 



inferred from the 24 — 70 fim colour assuming blackbody 
grains at a single temperature, Eq. 5). 

For extended disks r24-7o increases as r st i r moves 
outward. The disk T24-70 is located between the inner 
edge and where the surface density peaks, so moves out- 
ward with time. It is not simply at the peak surface den- 
sity because the inner disk contribution to the SED is 
non-negligible when r e a oc r 7//3 . Thus, for extended disks 
the increase in r24-7o with time is in fact slower than ex- 
pected by Equation (11). The fractional excess decreases 
somewhat because the decay of emission interior to r st ir is 
stronger than the increased emission from newly stirred 
regions. 

In contrast, belts show a small change in radius and 
then decay at near constant radius once stirred to the 
outer disk edge. The narrower disk extent means that 
7_4-70 is nearer the peak excess. When stirring reaches 
the outer edge, all disks break to a faster decline in / 
because no new regions can be stirred. Because the peak 
surface density is still evolving outward (i.e. these disks 
are evolving in the slow mode), T24-70 continues to in- 
crease. When t 3> t c at the outer edge, / decreases at con- 
stant T24-70 (because the disk has r e « oc r 7//3 everywhere 
and r24-7o cannot change). Because the collision time 
scales more strongly with radius than stirring (t c oc r 23 ^ 6 
vs. i s tir oc r 3 ), most disks do not reach evolution at con- 
stant T24-70 until long after the outer disk edge is stirred 
(and long after the maximum 1 Gyr age shown in Figure 
3). That is, slow evolution is more likely to occur at large 
radii (Eq. 10). 

For comparison, pre-stirred disks show similar trends 
to Figure 3, but / can only decrease because pre-stirred 
disks start evolving at t = everywhere. Thus, pre-stirred 
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disks do not show the initial increase in / seen for the 
self-stirred evolution. Also, there is no break to a faster 
decline in /, as occurs when self-stirring reaches the outer 
disk edge. 



4 APPLICATION TO A-STAR STATISTICS 

In this section we apply our model to evolution of 24 
and 70 zjm excesses around A-stars. We first outline the 
observations and the limitations of the pre-stirred model. 
Our aim is to test whether self-stirring can successfully 
reproduce the rise in 24 (im emission at ~ 10-30 Myr (and 
the 24 and 70 fim statistics at other times), and see what 
constraints the statistics set on model parameters. 

4.1 The A-star sample 

The observations are compiled from several sources. 
Ricke ct al. (2005) and Su et al. (2006) observed large 
unbiased samples of A-stars at 24 and 70 fim with 
Spitzer. While these samples provided the basis for the 
Wyatt et al. (2007b) study of pre-stirred evolution, they 
comprise stars mostly older than 10 Myr. 

To supplement these data, we collect 24 z^m data for 
B8-A9 stars in the following young clusters/associations: 
a Ori (Hernandez et al. 2007), OBla/b (Hernandez et al. 
2006), A Ori (Hernandez ct al. 2009), 7 Velorum 
(Hernandez et al. 2008), Upper Sco (Carpenter et al. 
2009), j3 Pic Moving Group (BPMG Rebull et al. 
2008), 6 NGC 2232 (Currie et al. 2008b), and IC 2391 
(Siegler et al. 2007). We use the K s -[24] excess to de- 
rive J^disk/^M,* for these objects and use spectral types 
from Kharchenko (2001) where needed. These data make 
a sample of about 400 A-stars observed at 24 (im with 
ages between ~3-800Myr. All stellar photospheres in 
these samples are detectable and there are no upper lim- 
its (i.e. disk detections are calibration limited). We do 
not include stars in more distant regions, where observa- 
tions are sensitivity limited and excess fractions are lower 
limits. 

Because our 24 /j,m sample includes stars younger 
than 10 Myr old, it is possible that some are proto- 
planetary disks. However, consistently shorter disk life- 
times around intermediate mass stars mean that few 
A-star protoplanetary disks survive beyond a few Myr 
(Kennedy & Kenyon 2009). Indeed we exclude only three 
stars; V346 Ori in OBla, HD 290543 in OBlb, and HD 
245185 in A Ori, which have K s -[24] « 6.5-7.5 and strong 
IRAC excesses (Hernandez et al. 2006). Less certain is 
whether the youngest stars in our sample contain dust 
left over from the protoplanetary disk phase, rather than 
from that created by fragmentation from self-stirring (e.g. 
Wyatt 2008). 

Figure 4 (left panel) shows the 24 fim excesses as 
a function of time binned into "large" (f24,tot/^24,* > 



2), "medium" (1.25 < fk.tot/Jk,* < 2), and "small" 
(f24,tot/-F24,* < 1.25) excesses. The age bins are 0- 
6 Myr, 6-20 Myr, 20-40 Myr, 40-190 Myr, 190-400 Myr, 
and 400-800 Myr. There are 106, 45, 37, 130, 51, and 
33 stars in these bins respectively. The primary result at 
24 fj,m is a decline in excesses on a ~ 150 Myr timescale 
(Rieke et al. 2005). 7 The large and medium excess bins 
add to give a peak in debris disk fraction at 30 Myr 
(e.g. Currie et al. 2008b). By separating the excesses into 
three bins, there is perhaps evidence for an additional 
trend: the large excesses peak around 10 Myr, and the 
medium excesses peak around 30 Myr. This behaviour 
suggests a population of large excess disks that form 
around 10 Myr, and then decay relatively rapidly to pro- 
duce the medium excess population a few tens of Myr 
later. 

Though these features are tantalising evidence of 
systematic trends that may be caused by self-stirring, 
their existence is motivated more by expectations than 
by statistical significance. In addition, Carpenter et al. 
(2009) note that because A-stars take ~10Myr to reach 
the main-sequence, stars of fixed spectral type decrease 
in luminosity from 1-10 Myr. The decreasing luminosity 
decreases D^i, and can cause at least some of the rise in 
24 /an excesses that has been attributed to self-stirring. 
Using the cumulative distributions shown in Figure 4 
(right panel) we estimate the significance of changes be- 
tween different age bins with the KS test. According to 
this test, each pair of cumulative distributions for the 
three youngest bins (three upper curves) could have all 
come from the same distribution (max KS probability of 
difference is 82%). There is at least a 99% chance that all 
other combinations of distributions are not drawn from 
the same distribution, with the exception of the fourth 
age bin when compared to the fifth age bin (78%). That 
is, there is a robust decay in excesses from early to late 
times. For the individual clusters collected above, aside 
from the 50 Myr IC 2391, there is no pair different at 
more than the 98% level according to the KS test, though 
NGC 2232 is different from o Ori, Orion OBlb, and 7 Vel 
at 93-98%, and Upper Sco is different from the BPMG 
at 92%. Therefore, while there are interesting trends at 
early times, the only formally secure result is that 24 /xm 
excesses decrease after ~50 Myr. 

While there are far fewer observations of A-stars 
at 70 fim, we use the Su et al. (2006) stars, comprising 
153 stars (including 19 observed with IRAS, and exclud- 
ing Be stars HD 21362 & HD 58715 and Herbig Ae/Be 
star HD 58647). These data include many upper limits, 
which are largely for stars farther than 100 pc. For the 
model comparison undertaken here, we use the sample 
as given by Table 4 of Su et al. (2006) (shown in Fig. 
9 of that paper, or Fig. 2 of Wyatt et al. (2007b)). As 
in Su et al. (2006) and Wyatt et al. (2007b), the data 
are binned into "large" (FVo,tot/FVo,+ > 20), "medium" 
(5 < F 7 o,tot/fYo,* < 20), and "small" (F 70 ,tot / 'F 70 ,+ < 5) 



6 Pic is removed from the Rieke et al. (2005) sample. We 
set the excesses of the Be stars HD 21362 (Rieke sample), HD 
67985 (7 Vol), and HIP 78207 (Upper Sco) to F 2 4,tot/-F 2 4,* = 
1. 



7 This decay is more obvious when viewed on a linear-time 
plot (e.g. Rieke et al. 2005; Su et al. 2006; Wyatt et al. 2007b) 
and is only seen in the oldest three bins in Figure 4. 
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Figure 4. Evolution of A-star 24 (im excesses. The left panel shows stars binned by age and large (triangles), medium (squares) 
and small (diamonds) excesses (ages offset slightly for clarity). Error bars are %/N in each bin and for the oldest large excess bin we 
assume 1 disk to calculate the error. Excesses reach an overall peak around 30Myr. The right panel shows cumulative distributions 
of 24 fim excesses for the same age bins, from dark to light grey lines. The dotted lines indicate the bin edges for small, medium, 
and large excesses. In the second youngest bin HR4796A and ft Pic are beyond the right edge of the plot. 




Figure 5. Same as Figure 4, but for 70 /jm excesses. 



excesses. The data are also binned by age: 0-30 Myr, 30- 
190 Myr, 190-400 Myr, and 400-800 Myr. 

The inclusion of upper limits in the 70 /im sample 
means that the excess fractions are upper limits. As a 
check, we compare this sample to an unbiased subset of 
stars closer than 100 pc (and exclude IRAS sources, all 
but 2/19 of which are upper limits). In this subset only 
two stars (HD 27962 & HD 142703) with upper limits fall 
above the small excess bin (with i*Vo,tot/-FVo,* < 5.61 & < 
7.63 respectively). The only significant difference between 
these two samples lies in the youngest age bin, where the 
<100pc subset has only one star (of 11) with a small 
excess, compared to 16 (of 31) for the full sample. This 
difference arises because the main sample includes many 
Upper-Sco sources (~150pc away) that have upper limits 
in the small excess bin. However, the paucity of young 
star forming regions within 100 pc also results in poor 
statistics in this bin (with only 11 stars). An alternative 
approach, setting all stars with upper limits to have small 
excesses, gives a ~10% increase in small excesses for the 



youngest two age bins, and little difference for the older 
two. For both approaches, there is little change in the 
excess fractions in all but the youngest age bin (relative to 
Poisson errors) . Based on these comparisons, we conclude 
that differences in the youngest age bin simply reflect 
poor statistics due to a lack of very close young stars. 
The small differences between (sub)samples in the older 
age bins mean that sample choice does not affect our 
results. Thus, we use the entire Su et al. (2006) 70 /mi 
sample for comparison with our model. 

Figure 5 shows the evolution of excesses at 70 /im. As 
with the older age bins at 24 /mi, these decay monotoni- 
cally with time, but with a slower timescale of >400 Myr 
(Su et al. 2006). A KS test shows that the differences be- 
tween adjacent age bins are not particularly strong (~90- 
98% chance of not being from the same distribution) , but 
it is very unlikely that non-adjacent bins could be taken 
from the same distribution (with probabilities >99.99%). 
Therefore, the decay from a higher to lower fraction of 
70 fim excesses over time is robust. 
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Finally, we use the Wyatt et al. (2007b) 46 star sub- 
sample of stars for which excess emission has been de- 
tected at both 24 and 70 fim by Spitzer (or 25 and 60 /an 
by IRAS). With the assumption of blackbody grains at 
a single temperature, this sample provides a set of disk 
radii (r24-7o) for comparison with our model. 

4.2 Method 

We now explore various population models and their abil- 
ity to reproduce the A-star statistics. For these models, 
we select 30,000 stars randomly in the spectral type range 
B8-A9 with uniform log(age), assuming all stars are born 
with planetesimal disks. As in Wyatt et al. (2007b) we 
fix some parameters and allow others to vary. We fix 
p = 2700 g cm -3 , q = 11/6, I = e and do not consider 
the effects of varying them. These assumptions leave , 
D c , e, rj, and the disk radii as parameters. 

In addition to comparison of fractional excesses in 
different age bins, our modelling process also involves 
comparison of model and observed disk radius distribu- 
tions Using T24-70- 

The low significance of the 24 /im trends means that 
we measure the ability of the self-stirring model to fit 
A-star statistics at early times partly with a formal x 2 
measure for the 24 and 70 /im statistics and the r24-7o 
distribution, but also the success of the model at repro- 
ducing trends in the data. In some cases we can rule out 
models based purely on a poor fit to the statistics for 
the older A-stars, in which case our conclusions are inde- 
pendent of the low significance of trends for the younger 
stars. 



4.3 Pre-stirred disks 

To illustrate the shortcomings of the pre-stirred model 
for young stars, the best fit parameters from Wyatt et al. 
(2007b) are compared to the A-star observations in Fig- 
ure 6. This Figure uses the model of §2 with t s ti r = 
so that disks are pre-stirred. The parameters are: = 
150 J kg -1 , D c = 60 km, e = 0.05, and a planetesimal 
belt radius distribution N(r) oc r 7 with 7 = —0.8 and 
disk width dr — r m id/2. We find that ?/ m id = 0.2 gives the 
best results (Wyatt et al. specified the disk mass distri- 
bution with M m id = 10 Mq). The Figure shows that the 
model still reproduces the 24 and 70 /im A-star statistics 
for stars 100 Myr and older, but fails to reproduce the 
trend in large and overall excesses for younger stars at 
24 /im. The model shows only a monotonic decline in the 
overall fraction of stars with excesses in contrast to the 
observations. 

As a metric for comparing models, we quote the x 2 
value at 24 and 70 /im, for the distribution of radii r24-7o, 
and the sum of all three. These are computed as the sum 
of squared differences between the data and model di- 
vided by the Poisson error for all 18 data points at 24 /im 
(six age bins and three excess bins), 12 points at 70 /im 
(four age and three excess bins), and 5 points for the 
?"24-70 distribution. In the pre-stirred case shown in Fig- 
ure 6, the values are X24 — 32.4, Xn = 5.2, and Xi = 17.0, 
for a total of Xtot — 54.6. Though these numbers give no 



indication of why a particular model succeeds or fails, or 
whether it reproduces the desired trends, they formally 
show how well the statistics are reproduced and provide 
a benchmark for measuring the success of different self- 
stirring models. 



4.4 Delayed stirring: the simplest case 

By taking the peak in large excesses at 10 Myr as an indi- 
cation of the delay time, we can add the simplest possible 
model of delayed stirring to our model. Figure 7 shows a 
model with almost the same parameters as the pre-stirred 
model with the addition of a 10 Myr delay (la width 
0.1 dex) before the onset of stirring for all disks (here 
??mid = 0.25). After this delay all disks begin evolving, re- 
gardless of radius. This simple modification to the model 
reproduces the 24 /im excesses better than the pre-stirred 
model with some change at 70 /im, yielding X24 = 19.9, 
X?o = 12.5, and x? = 21.9 for a total of Xtot = 54.3. 
The three 24 /im excess bins now show the right general 
trends: the large excesses peak at 10 Myr, the medium 
excesses show a broad peak over 10-100 Myr, and the 
small excesses reach a minimum around 10-30 Myr. The 
fraction of small 70 /im excesses at the earliest times in- 
creases due to more unstirred disks. An obvious physical 
motivation for this prescription is unclear, though there 
are some possibilities that may merit further study. 

One possibility might seem to be objects emerg- 
ing from a protoplanetary disk phase somewhat longer 
than the "typical" 6 Myr lifetime (e.g. Haisch et al. 2001). 
However, this explanation is unlikely for two reasons. Pri- 
marily, only a small fraction of stars have long lived disks, 
not the ~40% suggested by Figure 4. Further, primor- 
dial disk lifetimes around A-stars are consistently shorter 
than for less massive stars (Kennedy & Kenyon 2009). 
Therefore, the delay between primordial disk dispersal 
and the creation of large 24 /im excesses is too long for a 
plausible direct link. 

A possible delay mechanism is scattering of proto- 
planets and embryos to unstirred locations in the outer 
disk, where they stir smaller objects to collision veloci- 
ties. In this scenario the same fraction as in the 10 Myr 
peak, ~40%, of stars must harbour planets that undergo 
scattering events. This figure is higher than the frac- 
tions of intermediate mass stars both known (9% for 1.3- 
1.9 M Q , Johnson et al. 2007) and predicted (~20% for 
2-3 Mq, Kennedy & Kenyon 2008) to have gas giants. 
However, stellar systems may contain many more unde- 
tectable lower mass and/or more distant planets so the 
fraction of stars with planetary systems is probably not 
a major problem for this scenario. 

The main issues with a scattering scenario are tim- 
ing and depletion. For all scattering events to occur at 
around 10 Myr seems to produce a fine tuning problem, 
when simulations of scattering around Solar-type stars 
show that instabilities can happen at a wide range of 
epochs. For example, in the case of the "Nice" model 
for the Late Heavy Bombardment (LHB) in the Solar 
System, Gomes et al. (2005) show that the LHB epoch 
varies strongly with the primordial Kuiper Belt's inner 
edge location. In addition, this type of simulation is very 
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Figure 6. Best fit parameters for the Wyatt et al. (2007b) pre-stirred model at 24 and 70 (im (same symbols as Figure 4, joined 
by solid, dashed, and dotted lines for large, medium, and small excesses respectively). Symbols with error bars show observed 
A-star statistics. At 24 (im the model predicts too many stars with large excesses at the earliest times, and a monotonic decline 
in the overall fraction of stars with excesses. There are no new 70 fj,m data so the right panel is the same as Wyatt et al. (2007b) 
Figure 2. 
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Figure 7. Simple delayed stirring model compared with 24 and 70 fim statistics. Symbols and lines are as in Figure 6. All disks 
begin their evolution after ~ 10 Myr. 



sensitive to the initial conditions and it is unlikely that 
all systems undergo instability at any particular time. 

The way that planetesimal belts are depleted during 
scattering events also presents a problem. Based on an 
analysis of the Nice model, Booth et al. (2009) show that 
the signature of an LHB like event is a drop in 24 and 
70 /xm excess ratios of around 4 orders of magnitude with 
a brief (~ 15 Myr) peak in 24 /im excess, meaning that 
excess fractions should drop, not peak when LHBs occur. 
Therefore, planet-planet scattering is an unlikely source 
of the trends seen in the 24 A-star statistics. 

Another delay mechanism arises if planetesimals are 
initially large ("born big", Morbidelli et al. 2009). The 
lack of an initial population of small grains means that 
observable dust is not generated until the large planetes- 
imals collide and begin to decay, thus causing a rise and 
subsequent fall in the small dust population. This delay 
could plausibly be 10-30 Myr depending on the planetes- 
imal properties and locations. 



4.5 Self-stirring 

We now turn to the self-stirring model. As discussed in 
§3.1, D c is of order 10 km. The eccentricity when 1000 km 
objects form in the Kenyon & Bromley models is roughly 
constant with radius, and around e = 0.01. We therefore 
use these numbers as a order of magnitude guide and 
attempt to reproduce the A-star statistics with similar 
values. The remaining model parameters are Qn, r; m id, 
and the disk radii. For self-stirring, Zdciay has little im- 
pact on the results (see §3.2). Of particular importance 
is the inner hole specified by n n (or the minimum r m id 
for narrow belts), needed to ensure a peak in the overall 
fraction of stars with disks. 

We consider two types of disks in our population 
models. These disks are motivated by observations that 
show debris disks are commonly narrow belts, in contrast 
to the structure of protoplanetary disks that extend from 
several stellar radii to hundreds of AU. Though resolved 
observations may not discern between a truly narrow belt 
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and an annulus of bright emission in an extended disk, 
Figure 3 shows that the excesses for these disks evolve 
quite differently. Comparing population models to the A- 
star statistics may therefore suggest which kind of disk 
is more common. Because protoplanetary disks are the 
precursors of debris disks, knowing whether they have the 
same spatial extent will yield information about likely 
mechanisms and locations for planetesimal and planet 
formation. 



4-5.1 Narrow Belts 

Figure 8 shows the model that best reproduces the A-star 
trends for narrow belts. Belt radii are taken from a distri- 
bution with r m id = 15-120 AU and 7 = —0.8. Belts have 
width dr = r mid /2, = 45 J kg - , and ry mid = 0.15. 
We set D c — 10 km and e = 0.025, which are simi- 
lar to the values expected from the Kenyon & Bromley 
(2008) models. The fit has improved significance over the 
pre-stirred and simple delay models with X24 — 22.4, 
X?o = 15.0, X r = 2.1 and X L = 39.5. 

Initially we expected that the 10-30 Myr delay for 
the rise in excesses would lead directly to an inner disk ra- 
dius using Equation (9) for t Bt i T . Therefore, with 77 ~ 0.15 
we expect r m id. m in ~ 30 AU for stirring to occur around 
10 Myr (Eq. 9). However, setting r m id,mtn to values larger 
than sal 5 AU means that fewer disks are concentrated at 
small radii, and the distribution of inferred radii and the 
70 /im excesses become inconsistent with the observations 
(see §4.6). Because we cannot simultaneously reproduce 
the statistics and have a peak in large excesses at 10 Myr, 
it is actually the 70 /im excesses and r '24-70 distribution 
that sets r m id, m i n . The timing of the medium excess peak 
is partly set by Qy> because many of these disks decay 
from younger large excess disks. The relative fractions 
in the excess bins are set by the distribution of surface 
densities. 

Though a reasonable fit to the 24 and 70 /im statis- 
tics and overall distribution of radii is possible, there is 
some difference in the distribution of radii with time 
(Fig. 8, lower right). At 100 Myr, the lower envelope 
of model disk radii starts to increase from 15 AU to 
~60 AU by 800 Myr. This change is caused by the de- 
cay of smaller disks below detectable levels. Seven disks 
lie well away from the region populated by the model. 
The most discrepant disks are the same six noted by 
Wyatt et al. (2007b) when compared to the pre-stirred 
population model. All disks with radii less than 10 AU 
have ///max > 20 and may therefore be transient (HDs 
38678, 115892, 3003, and 172555). The potential influence 
of PR drag was offered as an explanation for the remain- 
ing two, HD 2262 and HD 106591, which we maintain 
here and also apply to HD 97633. 

As discussed in §3.2, we expect self-stirred debris 
disks to have stirring times similar to their collisional 
times (i.e. 1Z ~ 10 6 ). For the narrow belt model we find 
that the 1Z distribution peaks around 2 x 10 6 with nearly 
all disks within ±ldex. 

In summary, we find a reasonable fit to the A-star 
statistics for self-stirred planetesimal belts with physi- 
cally plausible parameters. The model does not show a 



strong peak in large 24 /im excesses at 10 Myr because 
the disk masses and minimum belt radii required mean 
that many of the smallest disks are stirred at their inner 
edges before they are observed. However, the model does 
reproduce all A-star statistics and trends at both 24 /im 
and 70 /im. 



4-5.2 Extended disks 

We consider three possible populations of extended disks, 
loosely motivated by potential processes that may shape 
protoplanetary and debris disks: (1) Disks with some 
fixed inner hole size of a few tens of AU that extend to 
some variable outer radius, perhaps due to their natural 
size or truncation by stellar encounters or companions, 
(2) disks with a fixed outer radius and a variable inner 
hole size, perhaps cleared by planets that form at vari- 
ous locations, and (3) disks with fixed inner and outer 
radii, similar in structure to young protoplanetary disks. 
We do not consider disks with both variable ri n and r out 
because these are effectively a mix of belts and extended 
disks and provide little useful information about which 
type of disk is more likely. We again use D c = 10km, 
and e = 0.01 as a starting point. When it is fixed, we set 
?"out = 150 AU. The variable inner or outer radii are cho- 
sen from a power-law distribution between the maximum 
and minimum with power law index 7. 

For the most extended disks, those with fixed ri n 
and r ou t, a reasonable fit to the 24 fim statistics can be 
found with ri n around 40 AU (X24 ~ 35), but the 70 /im 
statistics and radius distributions are significantly differ- 
ent (xw ~ 190 and Xr ~ 50). Similarly, for variable n n 
with 7 = —0.8 and fixed r out , good fits to the 24 /im 
statistics (X24 ~ 30) and radius distribution (xi ~ 5) 
can be achieved, but at the cost of a poor fit to the 
70 fxm statistics (xto ~ 120). The high x?o m both cases is 
due to ~40-50% of the model population having medium 
70 /im excesses at late (>300Myr) times. 

This high fraction can be understood by looking at 
how extended disks evolve, shown in Figures 1 and 3. 
While the 24 /im excess declines as stirring moves further 
out in the disk, 70 /im excesses for extended disks only 
show small decreases while the disk is still being stirred. 
The rj — 1/3 disk in Figure 1 (middle solid curves) has a 
24 /on excess >2 between 7-100 Myr, and would thus con- 
tribute to the large 24 /im excess peak in the population 
model. The 70 /im excess for the same disk rises above the 
medium excess ratio of 5 at 6 Myr, and remains higher un- 
til after 1 Gyr. Even when stirring reaches the outer disk 
edge, the strong radial dependence on the collisional time 
means that at ~100 AU the excess only decreases slowly. 
Therefore, extended disks with fixed r ou t that have large 
24 /an excesses at early times have medium (or large) 
70 /im excesses for long periods of time. This evolution 
means that the A-star statistics rule out self-stirred ex- 
tended disks with fixed outer radii that evolve like our 
model. This conclusion is independent of the significance 
of the 24 /im trends for the younger A-stars. 

Though fixed r ou t leads to problems with the 70 /im 
excesses at late times, this issue can be addressed by al- 
lowing r ou t to vary. With fixed r- m = 15 AU and the distri- 
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Figure 8. Best fit of the model population to 24 and 70 £im statistics with narrow belts. Top left and right panels: binned 24 (im 
and 70 fim evolution, Lower left panel: overall (dashed line) and detectable (dotted line) model T24— 70 (omitting the four possibly 
transient sources noted in the text), and Lower right panel: r24-7o vs. t compared to the 46 star subsample (open squares show 
three additional sources detected by Carpenter et al. 2009, HIP 76310, HIP 77911, and HIP 80088). 



bution of r out set with 7 = —0.8, the 70 (im statistics can 
be reasonably reproduced (xn ~ 6). However, this radius 
distribution has trouble reproducing the 24 /im statistics 
(X24 « 100) and the T24-70 distribution for the 46 star 
subsample, predicting too many disks with small radii 
(Xr ~ 13)- Attempts to remedy this problem by increas- 
ing the number of wider disks (increasing the power law 
index 7) results in essentially fixing r out , which leads to 
the previous problem of too many medium excess disks 
at 70 fim at late times. Thus, extended disks with fixed 
Tin also have trouble reproducing the A-star statistics. 

The conclusion that extended disks cannot repro- 
duce the A-star statistics may depend on the particular 
simplifying assumptions needed to make an analytic pop- 
ulation model. A possible key difference is illustrated by 
Figure 1: the Kenyon & Bromley (2008) models decay at 
70 pm more rapidly than ours and the excesses are ~2- 
5 times lower after 100 Myr, probably due to continued 
accretion and stirring (see §3.1). 

Whether continued accretion and stirring will lower 
the predicted 70 fim excesses for extended disks at late 
times without making the radius distribution signifi- 
cantly steeper (i.e. more disks at large r being unde- 
tectable) is not clear, but the comparison in §5.2.1 of 



Kenyon & Bromley (2009) gives some indication. Their 
model predictions of 70 /im excess fractions for disks 
with fixed n n and r out are several times higher than the 
Su et al. (2006) observations at the latest times, as we 
found in §4.5.2. Therefore, it appears that extended disk 
models that include continued accretion and stirring still 
evolve too slowly at 70 /J,m, providing further evidence for 
the belt-like nature of debris disks. 



4-5.3 Summary 

Our population model suggests that debris disks are more 
likely to be narrow belts than extended disks. The formal 
\ 2 for narrow belts is an improvement over the pre-stirred 
model, and reproduces all A-star statistics, including the 
rise and fall in 24 jj,m excesses. The inability of extended 
disks to reproduce the statistics is due to how they evolve. 
Disks with fixed n n predict too many disks at small radii 
and disks with fixed r out predict too many 70 /jm medium 
excess disks at late times. 

The reason belts work well appears to be because 
they restrict the evolution of excesses. The power-law 
distribution of r m id means that most disks are at rela- 
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tively small r, and the evolution is truncated when stir- 
ring reaches the outer edge (see Fig. 3). The T24-70 distri- 
bution is therefore a closer reflection of the initial power- 
law distribution, rather than being set by self-stirred disk 
evolution as for extended disks. 



4.6 Constraints on the best fitting model 

With the self-stirring model we have somewhat more 
power to constrain parameters than the Wyatt et al. 
(2007b) pre-stirred model. This ability arises be- 
cause we have estimates of D c and e based on the 
Kenyon & Bromley (2008) results. The remaining param- 
eters left to fit the statistics are Qq, 77 m id, and the disk 
radii parameters r m id, m in and 7 

For the belt model surface density, we find ?i m id = 
0.15 gives the best fit. Because we find that debris disks 
tend to be belts, our use of a surface density law has only 
a minor impact on our model. There is a factor 5 differ- 
ence in surface density from inner to outer disk edges, 
much smaller than the range of 77 we consider. Conse- 
quently, the model does not strongly constrain the initial 
surface density power-law index 5. For our best fit model 
we use S = 1.5, but find that populations with S — 1 
produce similar results. 

We find QJj = 45 J kg -1 gives the best fit to the 
peak in medium 24 /xm excesses, similar to that used 
for the comparison in §3.1. This value is best consid- 
ered an effective value for the evolution, since a range 
of are expected for different size planetesimals and 
at different radii. Our Qf) is reasonable for weak rock 
and ice for D ~ 10 km for the ~100-1000 ms _1 range 
of collision velocities (Leinhardt & Stewart 2009). The 
value of 45 J kg -1 is a third of that found for the best fit 
model in (Wyatt et al. 2007b), but as in that paper Qr> 
and e are degenerate. Thus, our results do not change if 
g*5/6 e -5/3 ^ moo L6hne et a| (20O8) show that the 

assumption of a single planetesimal strength is within an 
order of magnitude of a more complex model that in- 
cludes a size-dependent planetesimal strength (their Fig. 
11), with the largest differences occurring at late times 
>lGyr. 

Constraining the range of disk radii within the nar- 
row belt model is also possible to some degree. The model 
has trouble producing a stronger peak in large excesses at 
10 Myr while retaining a reasonable fit to the statistics. 
The model shows a slightly stronger peak in large excesses 
at 24 ^m when the minimum r m id is increased to 20 AU, 
but the fit to the 70 /im excess fractions becomes worse 
because the relative fraction of wider disks increases (for 
fixed 7). Therefore, at late times there are more disks 
with medium 70 (im excesses, essentially the same prob- 
lem faced for extended disks with fixed r out . We find it 
difficult to fit the A-star statistics with r m id,min > 20 AU 
(where X24 = 32.5, x?o = 24.8, Xt — 5.2 and total 
XL =62.5). 

On the other hand, decreasing the inner hole size 
actually betters the formal significance of the model fit 
to the A-stars. With r mid = 3-120 AU and 77 = 0.45 we 
find X 24 = 14.0, Xto = 6.6, = 10.6 and X tot = 31.2. 
This model shows a monotonic decline in excesses as 



for the pre-stirred model. However, in contrast to our 
favoured belt model, this model predicts a population of 
disks younger than 100 Myr with T24-70 less than 15 AU. 
This region is empty in the lower right panel of Figure 8 
(though we ignore disks we deem to be transient). Thus, 
while the typical minimum belt radius could lie between 
3-15 AU, we favour the model with minimum belt radii 
^15 AU. 

The power-law distribution of disk radii is fairly well 
constrained to 7 ~ —0.8. This constraint arises due to 
the belt-like nature of disks, which means that the model 
7 must be similar to the observed distribution. 

The chosen age distribution has a small effect on 
the large excess peak at 10 Myr. If we set the minimum 
age to 1 Myr (instead of 3 Myr based on a Ori stars) 
then the peak becomes slightly stronger. The difference 
arises because including younger stars results in more un- 
stirred disks (with small excesses) in the youngest age 
bin. Though this effect is minor, it shows that uncer- 
tainty in stellar ages can affect the results of the popula- 
tion model. 



5 DISCUSSION 

In the previous sections we have shown how self-stirred 
debris disks evolve. The model makes predictions, some 
of which can be compared to photometric observations, 
such as how disk radii inferred from blackbody models 
are distributed and evolve. 

The most basic prediction of this kind for pre- and 
self-stirred models is that the radial location of peak emis- 
sion should increase with time (Fig. 2). However, most 
surveys have failed to find any evidence for this trend 
(e.g. Su et al. 2006; Najita & Williams 2005). Rhee et al. 
(2007) find an apparent increase in disk radii inferred 
from IRAS colours. Unfortunately, all sources with disk 
radii greater than 100 AU — old stars largely responsible 
for the trend — have yet to be confirmed with new obser- 
vations (i.e. with Spitzer). This picture is also compli- 
cated by the expected increase in the lower envelope of 
radii at late times as close-in disks drop below detection 
limits (Fig. 8 lower right; Wyatt et al. 2007b). Therefore, 
our model shows that the expected increase in radii may 
not be as obvious as predicted by Equation (11), and that 
other predictions of self-stirring could be more useful. For 
example, another prediction of self-stirring is that disks 
with large radii for their age should have higher than 
average disk surface density, because these disks stir to 
large radii the fastest. 

However, it is important that these predicted trends 
are not just compared with photometric observations, be- 
cause disk models based solely on SEDs can be degener- 
ate and/or uncertain. Resolved imaging is necessary to 
confirm or correct SED-derived estimates. Imaging is also 
needed to test predictions such as the surface density pro- 
files shown in §3.2. Below we take a detailed look at a 
subset of resolved debris disks with the aim of compar- 
ing observed disk characteristics with those predicted by 
models. 

Also, we have not addressed an alternative possi- 
bility to self-stirring, that disks are stirred by secular 
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perturbations from planets not co-located with the disk 
(Mustill & Wyatt 2009). In this section we show that 
a model with t a tir set by planet-stirring can reproduce 
the A-star statistics and suggest that these planets could 
cause debris disks to be narrow belts. Because both mech- 
anisms can fit the A-star statistics, high resolution imag- 
ing is the best way to differentiate between self-stirring 
and planet-stirring. 

5.1 Self-stirring vs. planet stirring 

The fact that the A-star observations are well reproduced 
with a self-stirring model shows that this mechanism may 
be important for debris disks. However, self-stirring is not 
needed as an explanation for systems such as ft Pic where 
a planet is the likely stirrer. It is therefore important to 
predict features that allow the stirring mechanism to be 
identified. In the case of individual systems influenced by 
planets these features are well known, and stem largely 
from the same influence that causes planetesimal random 
velocities to increase and collisions to be destructive. Sec- 
ular perturbations both warp the disk (e.g. ft Pic) and 
cause it to be offset from the star (e.g. Fomalhaut). These 
features may be imaged directly, or an offset inferred from 
peri/apocenter glow (e.g. Wyatt et al. 1999). On shorter 
timescales, objects on unstable orbits too close to plan- 
ets are ejected, which can result in sharp disk edges (e.g. 
Fomalhaut). The remaining way of identifying whether 
a planet may be the stirrer is to detect it directly (e.g. 
Fomalhaut, and perhaps ft Pic). 

Though these features provide a way of inferring a 
stirring mechanism (and discovering planets), they can 
only be applied to individual systems. Distinguishing 
the dominant stirring mechanism at a population level 
is more difficult, because planet-stirring introduces yet 
more parameters to the model. To briefly look at whether 
a planet-stirred population model can reproduce the A- 
star statistics, we use the stirring time assuming internal 
perturbers (Mustill & Wyatt 2009) 



= 5 x 10" 



(1-40 

e p i 



3/2 



pl,Jup T pl 



(12) 



in Myr, where the 'pi' subscripts indicate planet prop- 
erties and r is the disk location where the stirring time 
applies. 

The increased number of parameters allows more 
flexibility in reproducing the observed A-star statistics. 
For example, the A-star statistics can be reproduced as 
well as in Figure 8 for narrow belts if we set the stirring 
time with planet properties: M p i = 0.5Mj up , e p i = 0.1, 
r-pi = r mid/3. That is, each belt is assumed to have 
a 0.5Mj up planet with eccentricity 0.1 located at one 
third of its average radius. 8 We assume that planets form 
early — during the protoplanetary disk phase — so the for- 
mation time can be ignored. For comparison with the pre- 
vious models, this planet stirred model has X24 = 19.7, 



8 Setting Equations (9) and (12) equal and solving for r p i 
suggests r p i <x v / r^id. However, the stronger scaling may be 
needed to account for the stirring time being independent of 
disk mass. 



X70 — 7.3, x? = 4.4 and Xtot = 31.3. This example is 
unlikely to be the only type of planet distribution that 
reproduces the observations, and shows that distinguish- 
ing between self-stirring and planet-stirring is not yet 
possible by this method. Future studies of this type can 
use distributions of known exoplanet properties as input, 
though these are only complete to ~5 AU, and disks may 
be perturbed by planets at much larger radii (e.g. Foma- 
lhaut). 

Some other interesting points can be made if disks 
are stirred by planets. The decreasing upper envelope of 
24 /xm excesses for A-stars suggests that the stirrers are 
not too far interior (or exterior) to the disk, because the 
strong radial dependence on the stirring time means that 
disks far from their planets will be unstirred early and 
then luminous at late times. This conclusion explains the 
success of the above example, where the planet location 
scales with the disk radius. 

Because the stirring time is set by planet properties 
and not the disk mass, disks can stir in the fast mode 
and no longer have a maximum surface density T e ff,max 
(see §3.2). Neighbouring planetesimal orbits can begin 
to cross at non-zero eccentricity as they precess, so the 
collision velocity steps from zero to the forced eccentric- 
ity times the Keplerian velocity when the disk is stirred 
(Mustill & Wyatt 2009). We can also derive a condition 
similar to (10), but now using Equation (12) for the stir- 
ring time. For an interior planet-stirred disk to evolve in 
the slow mode: 



Tin 



e p i 



(•-4) 
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This relation is qualitatively different to Condition 
(10), because the planet-stirring time increases more 
strongly with radius than the decay time (see Fig. 6 of 
Mustill & Wyatt 2009). Thus, disks are more likely to 
evolve in the fast mode at large radii, because the disk is 
stirred so late that it would have decayed earlier if it were 
pre-stirred. Parameters that shorten the planet-stirring 
time, such as higher e p i or M p \, or larger r p i (bringing 
the planet closer to the disk because this example is for 
interior planets) make the disk more likely to stir in the 
slow mode. 

In summary, a population model of disks stirred by 
secular perturbations from planets can reproduce the A- 
star statistics if the planets are located near the disk. 
This model is unlikely to be unique, as different distri- 
butions of planet properties can probably give similar re- 
sults. Therefore, it is not possible to distinguish between 
self-stirring and planet-stirring for debris disk popula- 
tions by this method yet. The features shown by high 
resolution imaging of individual objects, such as warps 
and offsets, remain the best marker of debris disks influ- 
enced by planets. Our planet-stirred example also moti- 
vates planetary system architecture as a possible reason 
for debris disks to be narrow belts. 
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5.2 The origin of narrow belts 

Returning to the idea that resolved disks can be roughly 
split into extended disks and belts (Kalas et al. 2006), 
our results suggest that debris disks tend to be nar- 
row belts that have minimum radii of ~15AU. A sim- 
ilar conclusion was reached by Chen et al. (2006), who 
found that most of their objects' IRS spectra were best 
fit by single temperature blackbodies colder than 130 K 
(thus also suggesting that disks have inner holes). In con- 
trast, nearly all young stars have evidence for protoplan- 
etary dust and gas disks that extend from very near the 
star (e.g. Haisch et al. 2001) to hundreds of AU (e.g. 
McCaughrean & O'dcll 1996; Watson et al. 2007). The 
implication is that not only the primordial disk extent 
sets where debris disks reside, but other influences such 
as photoevaporation, disk fragmentation, and truncation 
and clearing by planetary and stellar companions. Within 
the context of the previous section, planets that stir the 
disk may also be responsible for clearing it at other loca- 
tions. As is likely the case with /3 Pic (Augereau et al. 
2001), apparently extended disks may result from the 
blowout of small grains created in a relatively narrow 
planetesimal belt. 

To produce debris disks that are narrow belts, these 
mechanisms need to plausibly reproduce two qualitative 
trends: (1) disk inner and outer radii are positively cor- 
related because we find that they are narrow belts, and 
(2) most disks have relatively small radii, to reproduce 
the observed power-law distribution of disk radii. 

One possibility is that the belts are locations where 
planetesimal formation was possible or favoured. For ex- 
ample, one process that both clears dust from inner re- 
gions and enhances more distant regions is the influence 
of photoevaporative disk clearing. Alexander & Armitagc 
(2007) show that after the inner gas disk has cleared, 
small grains (<10-100cm) are dragged outward as the 
inner edge of the gas disk moves outward. At some point, 
the gas disk either becomes too tenuous to keep mov- 
ing the grains, or the dust density becomes comparable 
to the gas density. Either way, a concentrated mass of 
grains is left behind by the gas (though objects larger 
than ~lm are less affected by this process). Formation 
of Pluto-size objects will be enhanced here, either simply 
due to the faster growth time, or a more rapid instability 
(e.g. Youdin & Shu 2002). Another mechanism that may 
result in an inner hole and planetesimals at a particular 
location is the direct or rapid formation of planetesimals 
in the spiral arms of a self-gravitating disk (Rice et al. 
2006; Clarke & Lodato 2009). This process necessarily 
occurs beyond tens of AU where the disk is marginally 
stable and dust may be concentrated in spiral arms on 
a timescale shorter than the orbital period. This process 
may result in narrow belts, because at >100 AU distances 
the disk is unstable to fragmentation and may form com- 
panions that truncate the disk. However, given that only 
~10% of the Su et al. (2006) sample have known com- 
panions, binary truncation seems an unlikely process for 
setting disk outer radii. 

The alternative is that the belts are locations where 
systems are able to retain planetesimals. For example, 
one possibility is truncation by exterior stellar compan- 



ions, or within a cluster environment. However, stellar 
flybys are unlikely to be what sets disk outer radii be- 
cause the cross section for a close encounter suggests that 
this mechanism should more often result in large disks, 
rather than small ones. 

Planet formation provides a possible explanation of 
debris disk radii, and is naturally consistent with both 
self-stirring and planet-stirring models. It is reasonable to 
think that planetesimals form out to radii some fraction 
farther than where planets can form. Dynamical clear- 
ing by planets can then set disk inner radii, analogous 
to how Solar System planets dictate the Asteroid and 
Kuiper belt locations. In this picture debris disk systems 
therefore consist of an inner planetary system with some 
radial extent and a narrow planetesimal belt that extends 
somewhat further. This picture is essentially that of the 
planet-stirred example that reproduces the A-star statis- 
tics above in §5.1. 

The range and distribution of disk radii may be 
linked to the initial protoplanetary disk mass. Higher sur- 
face density disks are expected to form more giant planets 
over a wider range of radii (e.g. Kennedy & Kenyon 2008) 
and these systems are likely more susceptible to scatter- 
ing, resulting in more extended dynamical clearing and 
debris disks with larger inner radii (and provides a sim- 
plistic explanation for why the BPMG and TW Hydrae 
A-star disks with the largest excesses have the largest 
radii as discussed below). If most protoplanetary disks 
are relatively low mass (e.g. Andrews & Williams 2005), 
this scenario would also typically result in debris disks at 
relatively small radii in agreement with the T24-70 from 
our A-star sample. 

The degree to which planets influence debris disk 
structure probably varies. In some cases planets may sim- 
ply dynamically clear inner regions while the rest of the 
disk is self-stirred, whereas in other cases the disk struc- 
ture may be entirely set by migration and shepherding, 
scattering, and secular stirring. We now turn to a small 
sample of resolved A-stars that allow us to study these 
possibilities for individual systems. 



5.3 Comparison with resolved imaging 

Resolved debris disks that show structures such as warps 
and offsets reveal planets that may remain otherwise in- 
visible. In these cases, planet-stirring is probably more 
important than self-stirring. Therefore, resolved imaging 
allows estimation of the stirring mechanism in individual 
cases. 

Resolved imaging also allows other comparisons be- 
tween observations and models to be made, such as with 
the surface density profiles shown in Figure 2. Though 
the distribution of radii in our model is set by compar- 
ing model and observed r '24-70, this measure tends to 
underestimate disk radii and cannot account for disks 
with several dust components, providing further motiva- 
tion for imaging. 

As a sample of resolved debris disks, we use stars 
in the ~12Myr old f) Pic Moving Group (BPMG, 
Zuckerman & Song 2004) and the similarly aged ~8Myr 
old TW Hydrae Association. This sample is therefore 
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Table 1. BPMG and TW Hydrae A-stars with disks. Inferred 
disk radii T2A-io an< I / = ^disk/-^* from Wyatt et al. (2007b). 
Real disk radii r rea i are derived from imaging and modelling 
(J3 Pic and HR 4796A, Augereau et al. 2001; Telesco et al. 
2005; Schneider et al. 1999, 2009) or detailed SED modelling 
(Chen ct al. 2006; Smith et al. 2009). The fractional luminosi- 
ties i?24 = ^24,tot/p24,* are from Rcbull ct al. (2008) and 
Rieke et al. (2005). 



Name 


7*24-70 


''real 


^disk/^* 


#24 




(AU) 


(AU) 






HR 4796A 


27 


70 


330e-5 


97 


P Pic 


21 


70 


140c-5 


26 


HR 7012 


4 


0.9, 6 


50e-5 


5.9 


r) Tel 


25 


3.9, 24 


20e-5 


5.5 



roughly coeval and about the age of the 24 /im excess 
peak. 

The A-stars in this sample are HR 4796A (HD 
109573), p Pic (HD 39060), r] Tel (HD 181296), HR 
7012 (HD 172555), HR 6070 (HD 146624), and HR 
6749/HR 6750 (HD 165189/HD 165190). Of these, HR 
6070 and HR 6749/6750 have no excesses at 24 or 70 ^m 
(Rebull et al. 2008). Characteristics of the remaining four 
are shown in Table 1. The actual radii (r rea i) of HR 4796A 
and P Pic are several times larger than r24-7o, and in rea- 
sonable agreement for HR 7012 and the outer dust com- 
ponent of 7] Tel. The difference in the case of P Pic and 
HR 4796A is explained by the presence of small grains, 
which emit inefficiently at long wavelengths, and are thus 
hotter than a blackbody grain at the same stellocentric 
distance. 

Compared to the entire A-star sample, the 24 fim 
excesses in this 10 Myr sample are among the largest (see 
Fig. 4). Comparing the fraction of stars with disks with 
i?24 > 5.5 (i.e. the same or brighter than r\ Tel), the 
~ 10 Myr sample has 67 ±33% (4/6) while the fraction for 
the remaining stars in the 6-20 Myr age bin in the overall 
A-star sample is 8 ± 5% (3/38). Whether this relatively 
high excess fraction is the result of evolution within a low 
density association or simply due to a small sample size 
is unclear. 

Starting with P Pic, this star harbours the second 
brightest known debris disk, which appears extended 
in both scattered light and at IR and longer wave- 
lengths (e.g. Smith & Terrile 1984; Chini ct al. 1991; 
Holland et al. 1998; Telesco et al. 2005). The disk ap- 
pears over a wide range of radii, at least partly 
due to blowout of small grains created in collisions 
(Augereau et al. 2001). Modelling of the scattered light 
and mid-IR emission suggests that the parent disk where 
these grains originate is centred at around 70-80 AU 
(Augereau et al. 2001; Telesco et al. 2005). In addition, 
a ~10Mj up planet at ~10AU has been inferred as 
the cause of a warp in the disk (Mouillet et al. 1997) 
that reaches a peak near the same location as the 
parent belt (Heap et al. 2000). This proposed compan- 
ion has recently been imaged (Lagrange et al. 2009a), 
though confirmation with second epoch observations 
is needed (Lecavelier Des Etangs & Vidal-Madjar 2009; 
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Figure 9. Optical depth in the /3 Pictoris debris disk derived 
from Mid-IR imaging (Telesco et al. 2005). Ovcrplotted is the 
r c ff profile for a secular-stirred planetesimal belt (dashed line) 
between 40— 100 AU (see text for details). 



Lagrange et al. 2009b). Finally, Telesco et al. (2005) find 
an apparent clump of small grains in the South- West wing 
of the disk at ~50 AU. 

The solid line in Figure 9 shows the optical depth 
profile of p Pic derived from mid-IR imaging of P Pic 
(North-East wing, Telesco et al. 2005). The optical depth 
peaks at around 75 AU, with a decrease to larger and 
smaller radii. This profile suggests that the disk has been 
stirred to 75 AU. Comparison with Figure 2 suggests that 
the profile is more akin to a fast self-stirred disk because 
the profile drops immediately interior to 75 AU. That is, 
the collisional time is shorter than the stirring time at 
r st ir so Condition (10) is not satisfied. As discussed in 
§3.2 this profile is unlikely to be the result of self-stirring 
(though can be fit with a self-stirred model). 

There are, of course, other possible explanations for 
the origin of the p Pic optical depth profile. The sur- 
face density of small grains interior to 75 AU may drop 
more rapidly due to continued accretion and stirring by 
Pluto-size objects for example. Alternatively, if we still 
assume that the profile is due to collisional evolution, 
planet-stirring is a possible scenario. As noted is §5.1, this 
evolution can produce disks that stir in the fast mode. In- 
deed, a planet has already been proposed as the cause for 
a warp in the P Pic disk, and the maximum extent of the 
warp suggests that stirring due to the planet has reached 
~75 AU, the location of the peak optical depth. 

The dashed line in Figure 9 shows the optical depth 
profile of a planet-stirred model with good agreement 
between 40-80 AU. The model has a planet with mass 
M p i/Mq = 16 x 10~ 3 at dpi = 10 AU with e pl = 0.1, 
and the disk is therefore stirred to ~75AU in 12 Myr 
(though Eq. 12 shows that M p i, a p i, & e p i are degener- 
ate in setting £ st i r )- The value e = 0.025 is roughly the 
forced eccentricity set by the planet which in turn sets 
the collision velocities as ~150ms _1 . For the disk to be 
depleted interior to 75 AU we need fairly small planetes- 
imals, with D c = 0.1km, QJj = 40 Jkg -1 , and r\ = 0.01. 
To match the overall level of optical depth requires de- 
creasing q to 1.8 (from 1.83). The collision velocities may 
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be higher due to increased planetesimal inclinations, ex- 
pected if the planet is inclined relative to the disk, which 
is the interpretation of the observed warp (Mouillet et al. 
1997; Augereau et al. 2001). In this case, a model with 
larger planetesimals can reproduce the observed optical 
depth profile. 

To match the optical depth exterior to 80 AU re- 
quires very large x delay ~ 0.5. It is more likely that the 
emission outside 80 AU is due to small grains created in 
collisions at <80 AU being blown out of the system, thus 
making the disk appear more extended than it really is 
(Augereau et al. 2001). Because we infer small D c , an 
alternative explanation for the clump of small grains in 
the South- West wing at about 50 AU is needed, because 
the clump has about as much mass as a 100 km plan- 
etesimal (Tclesco et al. 2005). Collective phenomena such 
as resonance trapping (Wyatt 2003) or dust avalanches 
(Grigorieva et al. 2007) would be required to explain the 
clump. 

The disk around HR 4796A has a similar radius to 
/3 Pic. Detailed modelling suggests a ~15AU wide par- 
ent belt at 70 AU with a wider distribution of smaller 
blowout grains (Wahhaj et al. 2005). As with /3 Pic, this 
disk is unlikely to be self-stirred. The relatively sharp 
inner disk edge and a brightness asymmetry and pos- 
sible offset (e.g. Telesco et al. 2000; Wahhaj ct al. 2005; 
Schneider ct al. 2009), make HR 4796A reminiscent of 
the Fomalhaut disk, which is known to harbour an inte- 
rior planet that may affect the disk (Chiang et al. 2009). 
If the inner edge of the HR 4796A disk is truncated by a 
planet then it is probably stirred by that planet. 

The two remaining disks, those around HR 7012 
and rj Tel are fainter than the previous two, but still 
have large excesses relative to the overall A-star sam- 
ple. HR 7012 has a very small disk, with detailed mod- 
els of IRS spectra suggesting grain temperatures corre- 
sponding to 0.9-6 AU (Chen et al. 2006). These models 
require sub-/xm-sized grains, with a composition indica- 
tive of dust produced in a recent collision (Chen et al. 
2006; Lisse et al. 2009). The mass of grains inferred is 
of order 10 21 g, the mass contained in a planetesimal a 
few tens of km in diameter. This disk has an unusually 
small radius for its age (Wyatt et al. 2007b) and a rea- 
sonably high ///max ~ 100, also suggesting that the dust 
is likely transient and not due to self-stirring. However, 
this conclusion does not mean that the disk was never 
sclf-stirred. To have reached what may be an analogous 
stage to the giant-impact period that formed the Solar- 
System's terrestrial planets, objects orbiting HR 7012 al- 
most certainly went through the stages of growth where 
sclf-stirring is expected. This late stage of chaotic growth 
could be considered a second phase of self-stirring, where 
the stirring this time arises because the surface density of 
smaller objects is insufficient to damp the largest objects 
and big objects stir each other. 

Finally, r\ Tel has distinct planetesimal belts at 
~4 and 24 AU, each contributing about equally to the 
excess at 24/xm. (Chen et al. 2006; Smith et al. 2009). 
Smith et al. (2009) find that the outer disk can be ex- 
plained by a self-stirred model. If stirring has reached 
24 AU, the hotter dust may be the result of a collision as 
planets continue to grow in inner regions. In contrast to 



HR 7012 however, there is no evidence for sub-pm grains 
in the rj Tel disk from which it might be inferred that 
a recent collision is required to explain short lived dust 
(Chen et al. 2006). 

In summary, when confronted with detailed observa- 
tions, the self-stirring model appears to face competition 
from the continued growth of planets through stochastic 
collisions and their dynamical effects in trying to pro- 
vide explanations of disk structure. The /3 Pic disk is 
likely stirred by the proposed planet at ~10 AU. The non- 
azimuthal symmetry of the HR 4796A disk also suggests 
a planetary influence. The rj Tel disk is consistent with a 
self-stirring model and HR 7012 is probably a transient 
disk resulting from a recent collision. 



6 CONCLUSIONS 

Recent observations of young A-stars show evidence for 
an increase in the level and frequency of 24 jim. ex- 
cesses from ~3 to 10-30 Myr (Hernandez et al. 2006; 
Currie et al. 2008a, b). Excesses then decline on an 
timescale of ~150yr (Rieke et al. 2005). The rise in de- 
bris disk emission at early times has been interpreted as 
evidence for self-stirring, where a collisional cascade be- 
gins when Pluto-size objects form and stir planetesimals. 
Because the time taken to form Plutos increases with ra- 
dial distance from the central star, Wyatt (2008) noted 
that the 10-30 Myr delay also implied that A-star disks 
must have inner holes of order 10 AU if they are self- 
stirred. Though there is in fact little evidence that the 
fraction of stars with 24 jim excesses changes in the first 
50 Myr or so (Fig. 4) , the overall trend shown by the A- 
star statistics provides tentative observational evidence 
of self-stirring. However, a promising alternative to self- 
stirring should also be considered, that debris disks are 
instead stirred by secular perturbations from an eccentric 
planet (Mustill & Wyatt 2009). 

In this paper, we use the analytic model de- 
scribed in §2 to study the evolution of self-stirred 
disks. Our model is essentially the steady-state evolu- 
tion model of Wyatt et al. (2007b), modified to include 
sclf-stirring. We first compare our model to the detailed 
Kenyon & Bromley (2008) results, using an empirical de- 
lay for self-stirring to ensure we reproduce their excess 
trends over a range of disk surface densities and stellar 
masses (§3). The only difference in evolution is after the 
peak excess, with the Kenyon & Bromley models decay- 
ing more rapidly at 70 /im, probably due to continued 
accretion. 

We illustrate the implications of collisional evolution 
for resolved disks (§3.2). Because disks process their mass 
from the inside out, the surface density profile of any col- 
lisionally evolved disk region increases as r 7 ^ 3 and disks 
appear to increase in radius over time (pre- and planet- 
stirred disks show the same behaviour). The primordial 
surface density profile remains where the disk has not 
been stirred (and protoplanetary growth is ongoing). 

Disks with delayed stirring can evolve in two dif- 
ferent ways. If the collisional time is short compared to 
the stirring time (i.e. has more mass than it would if it 
were pre-stirred) then the disk rapidly loses mass as it 
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reverts to its equilibrium state. This evolution results in 
a bright narrow ring of emission where stirring is occur- 
ring (Fig. 2). While we suggest that this evolution is un- 
likely for self-stirred disks, it can occur for planet-stirred 
disks. Disks that stir before the collisional time evolve 
in the same way as pre-stirred disks, with the difference 
that there is less emission in exterior regions where the 
collisional cascade has not started. This is the typical 
evolution we expect for self-stirred disks. 

In §4, we turn to the observations and show why the 
pre-stirred debris disk model fails to produce the trends 
in the A-star statistics; with no mechanism to delay the 
onset of stirring, 24 /j,m excesses in the model are highest 
at the earliest times. The overall fraction of stars with 
disks declines monotonically with time in contrast with 
the observations, which peak around 30 Myr. Using the 
same power-law planetesimal belt radius distribution as 
Wyatt et al. (2007b), and planetesimal sizes and eccen- 
tricities consistent with the Kenyon & Bromley (2008) 
models, we show that the A-star trends and statistics 
can be reasonably reproduced by a self-stirring model 
with Qq = 45 J kg -1 and the average disk mass 0.15 
times an MMSN disk. Disks are "narrow belts" with 
width dr — r m id/2. The smallest planetesimal belt has 
r-mid = 15 AU, and the largest 120 AU. 

We have less success fitting the A-star observations 
with extended disks — disks with fixed inner and/or outer 
radii. Although the 24 /im emission can be reasonably 
reproduced with disks with fixed ~150AU outer radii 
and fixed or variable inner radii, these models result in 
too many 70 /xm excess disks at late times. This prob- 
lem arises because extended disks evolve at near con- 
stant 70 /im fractional luminosity until their outer edges 
are stirred (Figs. 1 & 3). Disks with fixed inner radii 
of ~15AU and variable outer radii also fail to fit the 
observed statistics, because the models over-predict the 
number of disks with small radii. Thus, our conclusion 
that debris disks are narrow belts and not extended is 
independent of the A-star trends for ages <50 Myr. 

Progress can be made in several directions to further 
understand the effects of self-stirring on model popula- 
tions. Our model only removes dust from the small end 
of the size distribution, whereas the Kenyon & Bromley 
(2008) models show that mass is also lost as Pluto-size 
objects continue to accrete fragments, and that mass loss 
is accelerated as the largest objects continue to grow. 
A model including continued accretion and stirring will 
lower 70 ^m excesses for the oldest disks, perhaps al- 
lowing extended disks to reproduce the A-star statistics. 
However, the Kenyon & Bromley (2009) model compari- 
son with A-star data suggests there will still be difficul- 
ties, with extended disk models that include continued 
accretion and stirring also evolving too slowly at 70 /im. 

Planets probably stir and set the structure of some 
disks. In §5.1 we show the A-star statistics can be fit with 
a population of narrow belts stirred by secular perturba- 
tions from an eccentric planet. The planet-stirred model 
produces essentially the same results as the self-stirred 
one, with the key to reproducing the A-star observations 
apparently being the planet location. If planets are lo- 
cated too far from the disk then stirring occurs too late 
and the characteristic ~150 Myr timescale decay of 24 ^m 



excesses does not occur. Thus, the successful model has 
0.5Mj up planets with e = 0.1 that are located at one 
third the disk radius. 

Therefore, population models cannot yet distinguish 
whether self-stirring or planet-stirring is more important. 
Population models are also unlikely to rule out one stir- 
ring mechanism due to the many model parameters. The 
poor statistics for A-stars younger than ~50 Myr also hin- 
der progress. The fact that the rise in 24 fim excesses for 
young A-stars has marginal statistical significance is un- 
likely to change in the near future, as most nearby regions 
have been studied with Spitzer and a significant increase 
in numbers awaits the launch of JWST. 

In §5.2 we consider the origin of narrow planetesimal 
belts, and suggest that planet formation provides a natu- 
ral explanation, if planetesimals form to radii somewhat 
larger than planets. The debris disk inner holes are then 
regions cleared by planets, and the outer extent set by 
where planetesimals can form. Depending on planetary 
system architecture, these planets may also stir the disk 
as suggested by our planet-stirred example in §5.1. 

In §5.3, we look more closely at the sample of 
~ 10 Myr old resolved disks around A-stars from the 
BPMG and the TW Hydrae Association, and find that 
only rj Tel allows a reasonable explanation with a self- 
stirring model. The disks around /3 Pic and HR 4796A 
seem more likely to be affected by planets. It is possi- 
ble that the /3 Pic debris disk is stirred through secu- 
lar perturbations from the planet proposed to orbit at 
~10AU. The disk around HR 7012 appears transient, 
though probably went through a phase of self-stirring 
when it was younger. These observations suggest that 
the degree to which debris disks are influenced by plan- 
ets varies, and that the answer to the question of debris 
disk stirring lies with high resolution imaging. 
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